Systematics of parton fragmentation in e + -e and nuclear collisions 
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Parametrizations of fragmentation functions (FFs) from e + -e _ and p-p collisions are com- 
bined with a parton spectrum model in a pQCD folding integral to produce minimum-bias 
fragment distributions. A model of in-medium FF modification is included. Calculated frag- 
ment distributions are compared with hard components from p-p and Au-Au p t spectra. Data 
are well described by pQCD over a large kinematic region for a range of Au-Au centralities. 
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I. INTRODUCTION 

RHIC collisions are conventionally described 
in terms of hydrodynamic (hydro) evolution of a 
thermalized bulk medium and energy loss of ener- 
getic partons (hard probes) in that medium. Hy- 
dro should dominate pt spectra below 2 GeV/c, 
parton fragmentation above 5 GeV/c, and "quark 
coalescence" in the intermediate pt interval. 

However, recent analysis of spectrum and 
correlation structure has revealed minijet struc- 
tures in RHIC collisions Q, 0, i, H H @, 0] ■ Two- 
component analysis of p-p and Au-Au spectra re- 
veals a corresponding hard component (minimum- 
bias fragment distribution), suggesting that jet 
phenomena extend down to 0.1 GeV/c [1, M\. 
Minijets [10|] appear to dominate the transverse 
dynamics of nuclear collisions at energies above 
\/snn ~ 15 GeV and provide unbiased access to 
fragment distribution structure down to a small 
cutoff energy for scattered partons (3 GeV) and to 
the smallest detectable fragment momenta (~ 0.1 
GeV/c). 

Minijets can be studied in the form of pt- 
spectrum hard components isolated via the two- 
component spectrum model. Measured hard 
components are compared with pQCD fragment 
distributions (FDs). Parton spectrum parame- 
ters and modifications to fragmentation functions 
(FFs) in more-central Au-Au collisions are in- 
ferred [Tl[]. The goal is a comprehensive pQCD 
description of all nuclear collisions. 



II. TWO-COMPONENT MODEL 

The two-component (soft+hard) spectrum 
model was first obtained from a Taylor-series ex- 
pansion of p-p p t spectra on uncorrected event 
multiplicity n c h for ten multiplicity classes 0]. 
The soft component was interpreted as longitudi- 
nal nucleon fragmentation, the hard component 
as transverse scattered-parton fragmentation. 

The two-component model for p-p collisions 
with soft and hard multiplicities n s + = n c h is 
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Coefficient nh/n s scales as ah^, So(yt) is a Levy 
distribution on m t and H^iyt) is a Gaussian plus 
QCD power-law tail on transverse rapidity yt = 
ln{(m t +pt)/mo} [8]. To compare with A-A spec- 
tra we define S pp = (1/yt) dn s /dyt with reference 
model n s Sq and similarly for H pp <-> Hq. 

The corresponding two-component model 
for per-participant-pair A-A spectra is 
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S NN {y t ) + v H AA (y t ;v) (2) 
SNN{yt) + vr AA (y t \ v) H NN (y t ), 



where Snn (~ S PP ) is the soft component and 
H AA is the A-A hard component (with reference 
H NN ~ H pp ) [|, HJ. Ratio r AA = H AA /H NN is 
an alternative to nuclear modification factor R AA . 
Centrality measure v = 2n binary / ' n partic i pant esti- 
mates the mean nucleon path length. We are in- 
terested in the evolution of hard component H AA 
or ratio r AA with A-A centrality. 
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FIG. 1: First: Fragmentation functions (FFs) from e + -e _ collisions for three energies with /^-distribution 
parametrizations (solid curves), Second: Full e + -e~ FF parametrization on parton rapidity y max , Third: FFs 
from p-p collisions for several dijet energies, Fourth: Full p-p FF parameterization on parton rapidity. 



III. FRAGMENTATION FUNCTIONS 



difference suggests that FFs may not be universal. 



e + -e~ (e-e) fragmentation functions (FFs) 
have been parametrized accurately over the 
full kinematic region relevant to nuclear colli- 
sions fl2[ ]. Light-quark and gluon fragmenta- 
tion functions D xx (x,Q 2 ) <-> D xx (y,y max ) (xx 
= e-e, p-p, A-A) are described above energy 
scale Q = 2 Ej e t ~ 10 GeV by a two-parameter 
beta distribution f3(u;p,q) on normalized rapid- 
ity u [12]. Fragment rapidity for unidentified 
hadrons is y = ln[(E + p)/m n ], and parton ra- 
pidity y max = ]n(Q/m n ). Parameters (p,q) vary 
slowly and linearly with y m ax 

above Q = 10 GeV 
and can be extrapolated down to Q ~ 4 GeV. 
Fig. Q] (first panel) shows measured FFs for 
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three energy scales from HERA/LEP [1 
The curves are (3(p, q) parametrizations which de- 
scribe data over the entire fragment momentum 
range. Fig. [1] (second panel) shows the FF en- 
semble vs energy scale Q as a surface plot [l2T |. 

Figure Q] (third panel) shows FF data from 
p-p collisions at FNAL The dotted line 

represents the lower limit for e-e FFs. There is a 
significant systematic difference between p-p and 
e-e FFs. The CDF FFs also reveal suppression 
at larger parton energies relative to LEP e-e sys- 
tematics. Fig. [T] (fourth panel) is a surface plot 
of the p-p FF parametrization [ll[ — the e-e FF 
parametrization modified by cutoff factor 



9cut 

(y) = tanh{(y - y )/£ y } y > y , (3) 

with yo ~ Cy ~ 1-5 determined by the CDF FF 
data [15]. The cutoff represents real fragment and 
energy loss from p-p relative to e-e FFs. The 



IV. pQCD FRAGMENT DISTRIBUTIONS 

The parton pt spectrum from minimum-bias 
scattering into an n acceptance near projectile 
mid-rapidity can be parametrized as 
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= fcutiVmax) A ymax exp{-(n QC D - 2) y m ax}, 

with y m ax = ln(2pt/m 7r ). The cutoff factor 

fcutiUmax) = {t&nh[(y m ax ~ y C ut)/£,cut] + l}/2(5) 

represents the minimum parton momentum 
which can lead to detectable charged hadrons as 
neutral pairs. Parton spectrum and cutoff pa- 
rameters are determined by comparing FDs with 
p-p and Au-Au spectrum hard components. 

Fig. [2] (first panel) shows the parton spec- 
trum (solid curve) with cutoff ~ 3 GeV inferred 
from a p-p pt spectrum hard component [ll[ . The 
bold dotted curve is an ab-initio pQCD calcula- 
tion [16]. The spectrum integrates to 2.5 ± 0.6 
mb, consistent with pQCD theory [l7( ]. 

The pQCD folding (convolution) integral 
used to obtain fragment distributions is 
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where D(y,y max ) is the FF ensemble from some 
collision system (e-e, p-p, A-A, in-medium or in- 
vacuum) , and da dijet / dy max is the parton spec- 
trum [111 ] . Hadron spectrum hard component 
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FIG. 2: First: Parton spectra inferred from this analysis for p-p collisions (solid curve) and central Au-Au 
collisions (dash-dotted curve) compared to an ab-initio pQCD theory result (bold dotted curve [HI), Second: 
pQCD folding-integral argument for e + -e~ FFs, Third: e-e FD (solid curve) and p-p hard-component reference 
(dash-dotted curve), Fourth: Fragment distribution (solid curve) compared to p-p hard-component data (points). 
Dotted curves correspond to ±10% change in parton spectrum cutoff energy about 3 GeV. 



d 2 rih/dy drj represents the fragment yield from 
scattered parton pairs into rj acceptance drj. Ef- 
ficiency factor e ~ 0.5 includes the probability 
that the second jet also falls within 5r]. At] ~ 5 is 
the effective 47r r\ interval for scattered partons. 
&nsd (~ 36 mb for A /s/v r /v = 200 GeV) is the 
cross section for NSD p-p collisions. 

Fig. [2] (second panel) shows integrand 



D ee {y,yr, 



da 



of Eq. ([6]) with unmodified 



max) dyn 

FFs from e-e collisions and lower bound at y m in ~ 
0.35 (p t ~ 0.05 GeV/c) (dotted line). Fig. M 
(third panel) shows the corresponding FD (solid 
curve), the "correct" FD describing inclusive 
hadrons from partons produced by free parton 
scattering from p-p collisions. The dash-dotted 
curve is the hard-component model inferred from 
p-p spectrum data 0. The FD from e-e FFs lies 
well above the measured p-p hard component for 
hadron p < 2 GeV/c (y < 3.3), and the mode is 
shifted down to ~ 0.5 GeV/c. The "correct" e-e 
FD strongly disagrees with the hard component 
of the p-p pt spectrum. Nevertheless, the e-e FD 
is the proper reference for nuclear collisions 111 ]. 

Fig. [2] (fourth panel) shows FD Hnn-voc as 
the solid curve, with measured FFs from p-p col- 
lisions. The mode of the FD is ~ 1 GeV/c. The 
solid points are hard-component data from p-p 
collisions and the dash-dotted curve is p-p model 
function H pp [8] . The comparison determines par- 
ton spectrum parameters y cu t = 3.75 (E cut ~ 3 
GeV), Ay max and exponent hqcd = 7.5 and es- 
tablishes a quantitative relationship among par- 
ton spectrum, measured FFs and measured spec- 



trum hard components over all pt, not just a re- 
stricted interval above 2 GeV/c. 



V. PARTON "ENERGY LOSS" MODEL 

Fragmentation in A-A collisions requires a 
model of parton "energy loss" or medium modi- 
fication. We adopt a minimal model of FF modi- 
fication (Borghini- Wiedemann or BW) [13] . Fig- 
ure [3] (first panel) illustrates the BW model (cf. 
Fig. 1 of [ll], £ p = \a{p jet /p) = ln(2p jet /m 7T ) - 
ln(2p/m n ) ~ y max — y). In- vacuum e-e FFs for 
Q = 14 and 200 GeV from the beta parametriza- 
tion are shown as dashed and solid curves fl2|. 
We can simulate BW accurately by changing pa- 
rameter q in (3(u'ijhq) by Aq ~ 1 (dash-dotted and 
dotted curves) [111 ]. Small reductions at larger 
fragment momenta (smaller £ p ) are compensated 
by much larger increases at smaller momenta. 
The largest changes (central Au-Au) correspond 
to an inferred 25% leading-parton fractional "en- 
ergy loss." Fig. [3] (second panel) shows the mod- 
ified e-e FF ensemble with FF modes shifted to 
smaller fragment rapidities y. 

Figure (third panel) shows H ee _ me( i (solid 
curve), the FD obtained by inserting in- medium 
e-e FFs from the second panel into Eq. ([6]). The 
dotted curve is the H ee _ vac reference from in- 
vacuum e-e FFs. The mode of H ee _ me( i is ~ 0.3 
GeV/c. Fig. [3] (fourth panel) shows results for 
p-p FFs. Major differences between p-p and e-e 
FDs appear below p t ~ 2 GeV/c (y t ~ 3.3). 
Conventional comparisons with theory (e.g., data 
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FIG. 3: First: e + -e~ FFs for two energies unmodified (solid and dashed curves) and modified to emulate parton 
"energy loss" [l8| (dash-dotted and dotted curves), Second: e + -e~ FF ensemble modified according to [HI], Third: 
Medium-modified FD from e + -e~ FFs (solid curve) compared to in-vacuum e + -e~ FD (dotted curve) Fourth: 
Medium- modified FD from p-p FFs (solid curve) compared to in-vacuum N-N FD (dotted curve). 



vs NLO FDs) typically do not extend below 2 
GeV/c figf | . The large difference between the two 
collision systems below 2 GeV/c reveals that the 
small-pt region, conventionally assigned to hydro 
phenomena, may be essential for effective study 
of fragmentation evolution in A- A collisions. 



VI. FRAGMENTATION EVOLUTION 

Measured FFs are combined with a 
parametrized pQCD parton spectrum to produce 
calculated FD XX for comparison with measured 
spectrum hard components H xx . Figured] (first 
panel) shows spectrum hard components Haa 
(solid curves) for five centralities from 200 GeV 
Au-Au collisions 0. The hard components 
scale proportional to nonary, as expected for 
parton scattering and fragmentation (jets). The 
points are from 200 GeV NSD p-p collisions 0]. 
The dashed curve is -ffNN-vacj and the upper 
dotted curve is -ff e e-med with Aq = 1.15, which 
corresponds to the most-central Au-Au curve 
(0-12%). The parton spectrum cutoff for i? e e-med 
has been reduced from 3 GeV (y ma x = 3.75) 
to 2.7 GeV (y 

max — 3.65) to match the central 
Au-Au hard component near yt = 3. 

Jet-related spectrum structure can also be 
studied with ratios. The conventional spec- 
trum ratio at RHIC is Raa- Because it in- 
cludes the spectrum soft component Raa strongly 
suppresses fragment contributions at smaller y t . 
Hard-component evolution with centrality is bet- 
ter resolved by ratio taa = Haa/H^n '■ However, 
studies in Ref. llj reveal that the proper refer- 



ence for all systems is the in-vacuum FD from 
e-e FFs, not p-p FFs. We therefore define ratios 
r xx = FD xx ^yyy/FD ee _ vac with xx = ee, NN, 
AA and yyy = med or vac to be compared with 
equivalent spectrum hard components H xx - yyy . 

Figure H] (second panel) shows ratios rede- 
fined in terms of the ee-vac reference: H pp (p-p 
data - points), Haa (peripheral Au-Au data - 
solid curve [9J) and calculated -ff e e-med (dash- 
dotted curve) and H^N-vac (dashed curve) all 
divided by reference H ee - Vac . Strong suppression 
of p-p and peripheral Au-Au data apparent at 
smaller yt results from the cutoff of p-p FFs. 

Figure U] (third panel) shows measured 
Hw/H ee-vac Ior more-central Au-Au collisions 
(solid curves) above a transition point on cen- 
trality at v ~ 2.5, with partial restoration of 
the suppressed region at smaller yt and strong 
suppression at larger yt- The latter has been a 
major observation at RHIC (high-p t suppression, 
"jet quenching" [20]). Newly apparent is the ac- 
companying large increase in fragment yield below 
2 GeV/c, still strongly correlated with the par- 
ent parton 0. Chan ges in fragmentation depend 
strongly on centrality near the transition point. 
It is remarkable that the trend at 10 GeV/c cor- 
responds closely to the trend at 0.5 GeV/c. H pp , 
Haa and ratios based on the e-e in-vacuum ref- 
erence are well described by pQCD FDs from 0.3 
to 10 GeV/c Hi). 
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y t y t y t 

FIG. 4: First: Hard-component evolution in central Au-Au collisions vs centrality 0. Large increases at smaller 
yt accompany suppression at larger y t . Second: FD ratios relative to an ee- vacuum reference for Au-Au collisions 
below the sharp transition, Third: FD ratios above the sharp transition revealing major changes in FD structure, 



VII. CONCLUSIONS 

Hard components of pt spectra can be iden- 
tified with minimum-bias parton fragmentation 
in nuclear collisions. Minimum-bias fragment 
distributions (FDs) can be calculated by fold- 
ing a power-law parton energy spectrum with 
parametrized fragmentation functions (FFs) de- 
rived from e + -e~ and p-p collisions. Alterations 
to FFs due to parton "energy loss" or "medium 
modification" in Au-Au collisions are modeled 
by adjusting FF parametrizations consistent with 
rescaling QCD splitting functions. The reference 



for all nuclear collisions is the FD derived from 
in-vacuum e + -e~ FFs. Relative to that reference 
the hard component for p-p and peripheral Au- 
Au collisions is found to be strongly suppressed for 
smaller fragment momenta. At a specific point 
on centrality the Au-Au hard component tran- 
sitions to enhancement at smaller momenta and 
suppression at larger momenta, consistent with 
FDs derived from medium-modified e + -e~ FFs. 

I thank the organizers of ISMD 2009 for a 
delightful and informative conference. This work 
was supported in part by the Office of Science of 
the US DOE under grant DE-FG03-97ER41020. 



[1] R. J. Porter and T. A. Trainor (STAR Col- 
laboration), J. P hys. Conf. Ser. 27, 98 (2005), 
hep-ph/0506172. 

[2] R. J. Porter and T. A. Trainor (STAR Collabo- 
ration), PoS C FRNC2006, 004 (2006). 

[3] J. Adams et al. (STAR Collaboration), Phys. 
Rev. C 73, 064907 (2006). 

[4] Q. J. Liu, D. J. Prindle and T. A. Trainor, Phys. 
Lett. B 632, 197 (2006). 

[5] J. Adams et al. (STAR Collaboration), J. Phys. 
G: Nucl. Part. Phys. 32, L37 (2006). 

[6] J. Adams et al. (STAR Collaboration), J. Phys. 
G 33, 451 (2007). 

[7] M. Daugherity (STAR Collaboration), J. Phys. 
G 35, 104090 (2008). 

[8] J. Adams et al. (STAR Collaboration), Phys. 
Rev. D 74, 032006 (2006). 

[9] T. A. Trainor, Int. J. M od. Phys. E 17, 1499 
f2008 UarXiv:0710.4504l 
[10] X. N. Wang, Phys. Rev. D 46, R1900 (1992). 



[11] T. A. Trainor, P hys. Rev. C 80, 044901 (2009), 

larXiv:0901.3387l 
[12] T. A. Trainor and D. T. Kettler, Phys. Rev. D 

74, 034012 (2006). 
[13] W. Braunschweig et al. (TASSO Collaboration), 

Z. Phys. C 47, 187 (1990). 
[14] M. Z. Akrawy et al. (OPAL Collaboration) Phys. 

Lett. B, 247, 617 (1990). 
[15] K. Goulianos (CDF Collaboration), Proceedings 

of the "QCD and high energy hadronic inter- 
actions," XXXII Rencontres de Moriond, Les 

Arces, France, March 22-29, 1997, 
[16] F. Cooper, E. Mottola and G. C. Nayak, Phys. 

Lett. B 555, 181 (2003). 
[17] I. Sarcevic, S. D. Ellis and P. Carruthers, Phys. 

Rev. D 40, 1446 (1989). 
[18] N. Borghini and U. A. Wiedemann, 

|hep-ph/ 0506218. 
[19] A. Adare et al. (PHENIX Collaboration), Phys. 

Rev. D 76, 051106 (2007). 



5 



Thomas A. Trainor 



[20] C. Adler et al. (STAR Collaboration), Phys. Rev. 
Lett. 89, 202301 (2002). 



6 



